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FOREWORD 


AFOSR  Grant  80-0136,  "Transonic  Wind  Tunnel  Measurements  on  a  Supercritical 
Wing"  was  initiated  by  the  Structures  and  Dynamics  Division  (FIB)  of  the  Air 
Force  Wright  Aeronautical  Laboratories  (AFWAL).  The  objective  of  this  effort 
was  to  obtain  experimental  data  on  a  supercritical  wing  to  guide  improvements 
in  transonic  unsteady  aerodynamics.  This  effort  was  a  part  of  a  cooperative  program 
involving  Lockheed-Georgia  Air  Force,  NASA,  and  NLR.  This  cooperative  effort 
has  become  known  as  the  LANN  program. 

Under  contract  F3361 5-79-C-3030,  NLR  conducted  aerodynamic  analyses 
and  correlated  with  the  available  test  data.  The  principal  investigators  on  the 
grant  at  NLR  were  J.  J.  Hors  ten,  R.  G.  den  Boer,  and  R.  0.  Zwaan.  The  principal 
investigators  on  the  contract  were  A.  Steiginga,  R.  Houwink,  and  R.  J.  Zwaan. 

The  AFWAL  project  engineer  was  L.  J.  Huttsell  and  the  work  unit  was  24010236. 

This  report  presents  the  correlation  of  theoretical  and  experimental  unsteady 
airloads  on  the  supercritical  transport-type  wing  {LANN  Model).  Parameters  in  the 
correlation  are  Mach  number,  frequency,  mean  angle  of  attack,  and  oscillation 
ampl itude. 
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INTRODUCTION 


During  the  first  part  of  the  cooperative  program  of  GELAC,  AFWAL, 
NASA-Langley  and  NLR,  unsteady  transonic  pressure  measurements  were 
performed  on  a  semi -span  wind-tunnel  model  of  a  transport -type  super¬ 
critical  wing,  oscillating  in  pitch  (Ref.  1 ) .  The  model,  of  which  the 
geometry  is  similar  to  that  of  wing  A  designed  by  GELAC  (Ref.  2),  has 
become  known  as  the  LANN  model.  During  the  second  part,  correlations  of 
experimental  data  and  analytical  results  from  calculations  were  to  be 
made  by  both  GELAC  and  NLR.  Eight  test  cases  with  attached  flow  were 
selected,  covering  influences  of  Mach  number,  frequency  and  mean  angle 
of  attack. 

The  present  report  contains  a  description  of  the  correlation  made  by  NLR 
and  presents  the  results. 

The  computational  method  used  by  NLR  is  the  so-called  quasi -3-d  (Q3D) 
method,  which  combines  2-d  unsteady  transonic  small  perturbation  theory 
(the  LTRAN2-NLR  computer  program.  Ref.  3)  with  2-d  and  3-d  subsonic  Doublet 
Lattice  methods  according  to  the  quasi -3-d  concept  (Ref.  It).  This 
engineering-type  method  was  developed  in  particular  for  applications  in 
flutter  calculations  for  high  aspect  ratio  wings.  A  first  application  to 
a  supercritical  wing  was  reported  in  reference  5. 

In  this  report ,  a  summary  of  the  LANN  model  measurements  is  described 
in  section  2.  An  outline  of  the  Q3D  method  is  given  in  section  3.  The 
computations  and  results  are  discussed  in  sections  ^  and  5,  respectively. 


2  EXPERIMENTAL  DATA 


2 

The  experimental  data  were  obtained  in  the  NLR  2  x  1.6  m  transonic 
windtunnel  (HST).  The  model  was  designed  and  built  by  GELAC.  During  the 
test  it  was  mounted  at  the  wind-tunnel  side  wall  and  could  pitch  about 
an  axis  normal  to  the  wall  through  62  %  of  the  root  chord  (Fig.  l).  Mean 
steady  and  unsteady  pressure  measurements  were  carried  out  at  six 
streamwise  sections,  using  many  pressure  tubes  connected  with  scanning 
valves  and  a  limited  number  of  pressure  transducers  mounted  flush  with 
the  model  surface  (Fig.  2).  In  four  other  sections  accelerometers  were 
mounted  to  measure  the  amplitude  distribution  of  the  oscillating  wing. 

The  Mach  number  was  varied  between  0.62  and  0.95  at  three  angles  of  attack, 

^Figures  are  located  at  end  of  report. 
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covering  a  considerable  part  of  the  transonic  flight  regime  (Fig.  3). 

The  Reynolds  number  was  approximately  5  x  10^,  based  on  the  mean  aero¬ 
dynamic  chord.  A  transition  strip  was  fixed  at  0.05  c..._  downstream  of 

MAC 

the  nose.  The  reduced  frequency  was  varied  between  k  »  0  and  0.3 
(f  =  0  to  72  Hz),  the  amplitude  was  =  0.25  deg  in  most  cases.  As  the 
lowest  natural  frequency  of  about  30  Hz,  corresponding  to  fundamental 
wing  bending,  was  well  within  the  measuring  range,  the  vibration  modes 
were  strongly  dependent  on  the  frequency  (Fig.  It).  The  experimental  results 
are  available  in  reference  1  as  unsteady  pressure  distributions,  sectional 
lift  and  moment  coefficients,  overall  lift  and  moment  coefficients  and 
normalized  vibration  modes. 

In  view  of  the  correlation  in  section  5  of  measured  and  calculated 
lift  and  moment  coefficients,  some  remarks  are  made  here  on  the  accuracy 
of  the  experimental  values.  Due  to  a  number  of  inoperative  pressure  holes, 
the  integrated  airloads  may  not  provide  an  accurate  approximation  of  the 
real  airloads  in  various  sections.  In  each  section  one  or  two  pressure 
holes  are  inoperative  in  the  first  half  of  the  sectional  chord.  For  flow 
conditions  where  the  shock  wave  is  close  to  the  inoperative  pressure 
holes,  the  unsteady  pressure  peaks  have  not  been  measured  accurately. 
Further,  inoperative  pressure  holes  are  also  located  in  the  nose  region 
at  the  lower  side  of  sections  5  and  6.  All  this  missing  information  about 
the  unsteady  pressure  distributions  may  have  serious  repercussions  for 
both  lift  and  moment  coefficients  and  may  lead  to  unrealistic  effects  of 
the  vibration  amplitudes  on  these  airloads. 


3  THE  COMPUTATIONAL  METHOD 

The  computational  method  used  at  NLR,  the  so-called  quasi -3~d  (Q3D) 
method  (Ref.  4)  is  of  an  engineering  type  and  has  been  developed  in 
particular  for  use  in  transonic  flutter  calculations  (Ref,  5).  The  method 
is  assumed  to  be  applicable  for  high  aspect  ratio  wings  for  which  spanwise 
gradients  in  the  wing  loading  are  small  over  a  considerable  part  of  the 
wing  span.  The  method  essentially  provides  spanwise  distributions  of 
unsteady  lift  and  moment  in  which  both  the  effects  of  finite  span  and 
wing  thickness  are  taken  into  account.  The  Q3D  method  is  based  on  the 
combination  of  three  methods: 
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-  2-d  and  3-d  Doublet -Lattice  methods  to  compute  subsonic  2-d  and  3-d 
sectional  unsteady  airloads, 

-  a  2-d  transonic  method  to  compute  2-d  sectional  unsteady  airloads 
including  the  effects  of  the  mean  steady  flow  field  and  of  shock  wave 
motions.  At  NLR  the  commonly  used  method  is  the  LTRAH2-NLR  code  {Ref.  3) 
which  solves  the  low  frequency  transonic  small  perturbation  potential 
equation.  Recently  also  a  viscous  version  has  become  available  (LTRANV 
program.  Ref.  6)  which  combines  the  LTRAN2-NLR  code  with  the  lag- 
entrainment  method  of  Green  for  a  steady  turbulent  boundary  layer. 


Spanwise  distributions  of  unsteady  airloads  are  computed  as  follows 
(Fig.  5): 

a.  The  wing  is  divided  into  a  number  of  streamwise  strips, 

b.  The  3-d  Doublet-Lattice  method  is  used  to  compute  sectional  unsteady 
airloads  for  unit  pitch  and  heave  motions  of  each  strip.  The  results 
form  a  matrix  [A^l  of  subsonic  3-d  aerodynamic  influence  coefficients 
(AIC's); 


c. 


d. 


A  diagonal  matrix  FA^pJ  of  subsonic  2-d  AIC's  is  generated  using 
results  of  a  2-d  Doublet -Lattice  method  for  unit  pitch  and  heave 
motions  of  each  strip,  taking  into  account  the  spanwise  variation  of 
reduced  frequency  based  on  the  local  semi-chord.  The  equivalent  2-d 
conditions  are  derived  from  the  3-d  streamwise  conditions  using  the 
well-known  equivalence  rules  for  a  swept  wing.  The  calculated  unsteady 
airloads  are  transformed  back  to  quantities  based  on  the  3-d  stream- 
wise  conditions. 

Analogous  to  the  subsonic  2-d  computations,  a  similar  matrix  of 

transonic  2-d  AIC's  is  generated  using  results  of  a  2-d  transonic 
theory.  In  order  to  reduce  computational  cost,  transonic  flow 
computations  are  usually  carried  out  for  one  representative  wing 
section  (for  the  LANN  wing,  section  3  at  1*7. 5  %  of  the  span).  The 
initial  steady  flow  conditions  for  the  unsteady  transonic  flow 
computations  are  chosen  such  that  the  computed  sectional  mean  steady 
pressure  distribution  can  be  matched  satisfactorily  to  the  experimental 
distribution,  taking  due  account  of  the  shock  wave  location.  The 
approximation  made  by  neglecting  the  variation  of  the  mean  steady 
pressure  distributions  over  the  wing  span  is  considered  to  be  reasonable 
in  view  of  the  approximative  nature  of  the  various  steps  in  the  Q3D 
method. 
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e.  The  results  are  combined  into  a  matrix  of  AIC's  for  the  3-d  transonic 
unsteady  airloads  as  follows : 


(1) 


The  right-hand  side  of  Eq  ( 1 )  describes  in  a  condensed  notation  the 
following  correction  procedure  for  the  sectional  unsteady  lift  and 
moment  coefficients  due  to  unit  motion: 


,  T  .  ,  T  \  ,  S  „  S 
( m . .  +  b. .  k. .  ;  k .  .  /k . . 

«3J2d  1J  jj2D  1J3D  j,32D 


(2) 

(3) 


where  b  is  the  shift  of  the  aerodynamic  center  due  to  3-d  effects: 


b.  . 
ij 


S  t.  S  S  /.  s 

m.  .  /k .  .  -  m. .  /k .  . 

IJ3D  iJgp  jj2D  j,32D 


<10 


f.  Expressing  the  spanwise  variation  of  the  vibration  modes  in  a  matrix 

T 

[q] ,  the  distributions  of  lift  and  moment  coefficients  [C  ]  are  given 
by: 

K]  ■  [A3d]M‘  <51 

g.  Approximate  3-d  sectional  transonic  pressure  distributions  are  obtained 
by  multiplying  the  2-d  transonic  results  for  the  local  motion  by  the 
local  ratio  of  3-d  to  2-d  subsonic  unsteady  lift  coefficients. 


l|  COMPUTATIONS  FOR  THE  LANN  MODEL 


A  program  of  eight  computational  test  cases  was  agreed  with  GELAC. 
This  program  is  presented  in  table  1.  The  corresponding  experimental  data, 
together  with  definitions  and  sign  conventions  are  derived  from 
reference  1  and  included  in  Appendix  A.  The  mean  steady  state  conditions 
are  a  subsonic  case  and  three  transonic  cases.  An  impression  of  the  types 
of  flow  is  given  in  figure  6  by  the  steady  pressure  distributions  on  a 
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TABLE  1 

LANN  unsteady  test  cases  for  computations 


Run 

nr- 

M 

00 

a 

m 

(deg) 

cz 

f 

(Hz) 

k 

a1 

(deg) 

effect  of: 

M 

i 

x> 

at 

m 

k 

al 

22 

0.62 

0.6 

0.28 

24 

0.099 

0.25 

1 

f 

24 

0.62 

0.6 

0.28 

48 

0.199 

0.25 

J 

i 

118 

0.77 

0.6 

0.30 

24 

0.080 

0.25 

i 

12U 

0.77 

2.6 

0.48 

24 

0.080 

0.25 

l 

i 

f 

126 

0.77 

2.6 

0.48 

48 

0. 160 

0.25 

1 

1 

73 

0.82 

0.6 

0.32 

24 

0.076 

0.25 

r 

79 

0.82 

0.6 

0.  32 

24 

0.076 

0.5 

1 

85 

0.82 

0.6 

0.32 

48 

0. 151 

0.25 
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representative  wing  section  (no  3  at  1+7.5  %  semi-span}.  The  eight 
unsteady  flow  test  cases  were  selected  as  representative  for  the  measured 
effects  of  Mach  number,  frequency,  mean  angle  of  attack  and  amplitude. 
Because  the  theoretical  methods  investigated  are  applicable  to  attached 
flow,  no  conditions  with  separated  flow  were  selected. 

The  unsteady  airloads  were  computed  with  the  3-d  subsonic  Doublet - 
Lattice  method  as  well  as  the  transonic  Q3D  method.  Both  sets  of  results 
are  presented  in  this  report.  In  the  latter  method,  the  2-d  transonic 
unsteady  flow  computations  for  section  3  were  carried  out  using  the 
inviscid  LTRAN2-NLR  code  and,  for  one  case,  also  the  viscous  LTRANV  code. 
The  streamwise  coordinates  of  section  3  are  listed  in  table  2.  For  the 
equivalent  2-d  computations,  the  2-coordinates  were  smoothed  and  scaled 
by  the  equivalence  rule  for  a  swept  wing  ( ze=z/cos  A  ) .  The  transonic 
unsteady  pressure,  lift  and  moment  coefficients  (Cp.  ,  C  and  C  ,  defined 

1  Z .  JS . 

in  Appendix  A)  were  derived  from  the  first  harmonic  components  ot  the 
computed  results  of  the  second  cycle  of  a  sinusoidal  oscillation.  Higher 
harmonics  were  very  small,  except  near  shock  waves. 

The  unsteady  computations  were  complicated  by  the  fact,  mentioned  in 
section  2,  that  the  measured  vibration  mode  varies  with  frequency  and, 
though  less  strongly,  with  the  mean  steady  state  conditions  (Fig.  1+). 

Also,  each  measured  vibration  mode  had  out-of-phase  components.  In  the 
present  computations  both  effects  were  taken  into  account.  The  effect  of 
the  out -of -phase  component,  however,  appeared  to  be  small  in  all  cases 
considered. 

In  addition  to  the  unsteady  flow  computations,  also  quasi-steady 
airloads  were  computed  for  each  mean  steady  state.  A  comparison  with 
quasi -steady  experimental  data  gives  a  useful  supplement  to  the  unsteady 
results,  because  no  errors  associated  with  the  tube  transfer  functions 
and  the  measurement  of  the  vibration  mode  are  introduced. 
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TABLE  2 


Coordinates  of  the  LANN  wing  model  section  3 
(n  ■  0.475,  Local  chord  =  258.06  mm) 


upper 

side 

lover  side 

x/c 

z/c 

x/c 

z/c 

x/c  z/c 

x/c  z/c 

x/c 

z/c 

n  *00000 

o.n 1001 

0.02822 

0,03617 

n.oooon  o.oiooi 

0,04299  -0,0205ft 

0*91418 

-o.n  tons 

0*00001 

0.01043 

0,031 14 

0.03717 

0.00008  0.00369 

0,04542  -0,02110 

0,9141ft 

-0,00963 

0*00019 

0,01 308 

0,03310 

0.01784 

0.00018  0.00515 

0,0479ft  -0,02706 

0.94849 

-0.00874 

0*00053 

0,01281 

0.03566 

0.01868 

0.00048  0.00394 

0.05045  -0.02275 

0,96149 

-0.01006 

0.00089 

0,01466 

0,01851 

0.03957 

0.00088  0.00305 

0,05315  -0.02349 

0. 96965 

-0.010*2 

n  *001  33 

0.01589 

0.04086 

0.04025 

0.00100  0.00282 

O.057H5  -0,02477 

0,97750 

-0.01089 

mom  *6 

0,01658 

0*04147 

0,04095 

0.00141  0.00201 

0,06279  -0.02605 

0.98324 

-0.0 1 113 

0.00204 

O.OL  741 

0.04668 

0.0417ft 

0.00172  0.00142 

0.06792  -0.02736 

0.9B728 

-0.0  I  169 

0,00220 

O.OI784 

0.04886 

0,04227 

0.00210  0.00083 

0 , 0  7  2  ft  i  -0 , 0  2  ft  5  2 

0*99289 

-0.0122* 

0,00287 

0.0 1 R55 

0.05101 

0.0427ft 

O.OO255  0.00014 

0.07745  -0.02963 

0 , 995flft 

-0.01255 

0.003 ! 3 

0*01933 

0,05309 

0,0432! 

0.00278  -0.00002 

0,08323  -0,03098 

1 .00000 

-0.01 298 

0.00338 

0,01943 

0.05520 

0.04365 

0.00328  -0.00078 

0.08734  -0.03188 

0.00184 

n.020]  i 

0,05646 

0,04392 

0.00345  -0.00090 

0,09256  -0,03303 

0,00438 

0,02105 

0  .05  7  76 

0,04416 

0.00391  -0.00150 

0. 09715  -0,03402 

0.00478 

0.02150 

0.06279 

0,04512 

0.00429  -0.00192 

0. 10215  -0.03507 

0.00494 

0.02164 

0,06762 

0.04593 

0.00469  -0.00237 

0,10725  -0.03613 

0.00533 

0.02207 

0,07253 

0.04670 

0.00525  -0.00284 

0.1 | 690  -0,03807 

0.00563 

0,02248 

0,07749 

0.04739 

0.00560  -0.00326 

0.12674  -0.03999 

0,00619 

0.02312 

0.08732 

0.04861 

0.00615  -0.00386 

0,13650  -0 . 0  4 l ft  5 

0.00657 

0,02339 

0.09770 

0,04975 

0.00644  -0.00412 

0.14639  -0,04366 

0.0069ft 

0.03392 

0,10686 

0.05065 

0.00652  -0.00404 

0,15625  -0.04541 

O.Of172  3 

0,02407 

0.11705 

0,05152 

0.00716  -0.00483 

0,16622  -0,04709 

0.00766 

0.02469 

0,12658 

0.05226  1 

0.00770  -0.00531 

0. 1761  I  -0.04866 

0.00821 

0.02522 

0,15216 

0,05387 

0.00815  -0.00563 

0,18594  -0.05015 

O.OOB53 

0,025  5  2 

0.1765 1 

0.05499 

0.00861  -0.00600 

0,1964!  -0,05163 

0,0090! 

0,02593 

0.20061 

0.05585 

0.00907  -0.00635 

0 . 2 05 ft  4  -0 . 0 5  2 ft  9 

0,00953 

0,02638 

0.22548 

0.05640 

0.00961  -0.O0675 

0.21521  -0,05407 

0,00964 

0,02644 

0,25834 

0.05683 

0.01011  -0.00702 

0.22575  -0*05533 

0.01004 

0,02672 

0*11258 

0,05679 

0.01069  -0.00741 

0,25009  -0.057*4 

0.0  1068 

0.02730 

0,37718 

0,05573 

0.01  l  10  -0.00769 

0.27463  -0,06020 

n.oim 

0.02761 

0,42746 

0,05418 

0,01 153  -0.00794 

0 , 2  9  ft  ft  ft  -0,06197 

0,01  i  9  5 

0.02819 

0,47585 

0.0520! 

0.01 249  -0.00848 

0.32343  -0,06314 

O.OI25 J 

0,0285ft 

0. 521ft? 

0,04929 

0.01155  -0.00901 

0 , 3  4  ft  9 i  -0,06433 

n.m  324 

O.O3908 

0,57090 

0.04580 

0.01456  -0.00951 

0,17271  -0,06495 

0,01381 

0,02947 

0.62482 

0.04129 

0.01540  -0.00990 

0.3997ft  -0,06522 

0.01 442 

0,02980 

0.6720J 

0.03658 

0.01642  -0.01040 

0,42228  -0,06507 

0,01 492 

0,01007 

0,71828 

0,03132 

0.01739  -0.01088 

0.44655  -0.06457 

0.01546 

0.01040 

0.77026 

0,02474 

0.01841  -0.01 135 

0,47] 14  -0.06361 

0.01601 

0.03076 

0,81917 

0.01771 

0.01950  -0.01 187 

0,4961?  -0.06232 

0.01630 

0.01085 

0.86493 

0.01053 

0.02041  -0. 01210 

0.52041  -0.06061 

0.01687 

0,03119 

0,91586 

0.00224 

0,02144  -0.01282 

0.55624  -0,05725 

0.01 745 

0.03155 

0,92894 

0.00000 

0.02235  -0.01121 

0.58729  -0.05374 

0*01800 

0,03182 

0.96489 

-0,00594 

0.02 142  -0.01 168 

0,6201?  -0.04935 

0,01844 

0. 01204 

1  .00000 

-0.01  165 

0.02414  -0.01 405 

0.65092  -0,04474 

0 .01960 

0,03256 

0.02517  -0.01643 

0 , 6  ft  2 1  ft  -0,01987 

0.02036 

0.03292 

0.02659  -0.01487 

0,71732  -0,03435 

0.02129 

0,03341 

0.02732  -0.01517 

0 .750 ft  3  -0,02ft 9  8 

0.02227 

0,03385 

0,02843  -0.01 556 

0,76812  -0.02599 

0,02328 

0.03426 

0.0307  9  -0.01 64  1 

0,79] 17  -0.02259 

0,02428 

0,01469 

0.03318  -0.01719 

O.ftHft!  -0.01915 

0,02553 

0,01518 

0.01561  -0.01 809 

0-84046  -0,01603 

0.02637 

0.03552 

0.01810  -0.01904 

0 ,864  81  -0.01140 

0,07731 

0,03587 

0.04044  -0.01979 

0,8894  6  -0,0(136 
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PRESENTATION  AND  DISCUSSION  OF  RESULTS 


5*1  Steady  pressure  distributions 

In  figure  6  the  computed  mean  pressure  distributions  in  section  3 
are  compared  with  experimental  results.  To  obtain  a  reasonable  matching 
of  theoretical  and  experimental  results,  the  theoretical  angles  of  attack 
differ  from  the  experimental  values.  This  is  shown  in  table  3,  where  the 
equivalent  2-d  values  of  and  used  in  the  computations  are  compared 
with  the  experimental  3-d  values.  With  increasing  Mach  number,  the 
theoretical  incidence  has  to  be  decreased,  to  prevent  the  shock  wave  from 
moving  too  far  downstream.  In  addition  to  the  neglect  of  3~d  effects  and 
viscous  effects,  the  used  small  perturbation  theory  may  explain  this 
phenomenon. 

As  a  result  of  the  matching  procedure,  a  reasonable  agreement  is 
obtained  between  theoretical  and  experimental  results.  The  calculated 
values  at  the  rear  part  of  the  chord  are  overpredicted  due  to  the  neglect 
of  viscous  effects.  In  the  transonic  results  the  pressure  level  in  the 
supersonic  flow  region  and,  consequently t  the  strength  of  the  shock  wave 
are  underpredicted.  This  also  points  to  shortcomings  of  the  small 
perturbation  theory, 

5.2  Unsteady  airloads  as  function  of  Mach  number 
5.2,1  Quasi-steady  results 

In  figures  7  to  10  quasi-steady  results  are  presented  for  =  0,62, 
0,77  and  0,82  at  a =  0,6  deg.  The  corresponding  mean  steady  pressure 
distributions  are  shown  in  figure  6a,  b  and  c. 

Pressure  distributions :  Figure  7  presents  the  quasi-steady  pressure 
distributions.  The  theoretical  results  give  a  good  qualitative  prediction 
of  the  development  and  downstream  displacement  of  the  shock  pressure  peak, 
with  increasing  Mach  number.  In  the  subsonic  case  a  good  agreement  with 
experiment  exists  except  for  the  nose  region,  where  much  higher  pressure 
peaks  are  predicted.  With  increasing  Mach  number  the  pressure  level  on 
the  front  half  of  the  airfoil  is  increasingly  overpredicted.  This  may  be 
largely  due  to  the  small  perturbation  assumption  underlying  the  2-d 
transonic  theory,  as  well  as  to  the  neglect  of  viscous  effects. 

It  should  be  noted  that  the  experimental  results  at  ~  0,77  were 
obtained  for  an  amplitude  of  a  -  1  deg,  which  leads  to  a  lower  and  wider 
pressure  peak  than  for  =  0,25  deg.  This  effect  of  amplitude  is  very 


8 


TABLE  3 

Equivalent  2-d  flow  conditions  for  matched  theoretical  and 
experimental  pressure  distributions  on  wing  section  3 


Run 

nr. 

M 

CO 

a 

m 

(deg) 

M 

09 

('M^cos  A  ) 

a 

m 

e 

(des) 

22,  2k 

0.62 

0.6 

0.562 

0.5 

118 

0.77 

0.6 

0.698 

0 

\2k,  126 

0.77 

2.6 

0.698 

0.75 

73,  79,  85 

0.82 

0.6 

O.T^ 

-0.25 
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clearly  illustrated  for  =  0.82,  where  the  most  significant  changes  in 
pressure  {occurring  at  the  shock  peak)  are  shown. 

Sectional  coefficients:  The  spanwise  distributions  of  quasi-steady  lift 

and  moment  coefficients  are  shown  in  figures  8  and  9. 

In  the  subsonic  case  (Mm  =  0,62)  the  computed  effect  of  thickness 

on  the  lift  coefficient  is  small,  as  expected,  and  consists  of  a  small 

increase  of  C'.  The  level  of  Cl  is  overpredicted,  which  can  be  largely 
z  z 

attributed  to  the  neglect  of  viscous  effects.  The  sectional  moments  are 
predicted  very  well. 

With  increasing  Mach  number,  the  computed  steady  flow  effects  in  C'z 

also  increase.  This  can  be  associated  with  the  development  of  shock 

waves  and  the  increased  extent  of  the  supersonic  flow  region,  as  shown 

for  section  3  in  figure  6a,  b  and  c.  Although  the  experimental  results 

also  exhibit  a  stronger  effect  of  Mach  number  on  C'  than  predicted  by 

z 

subsonic  theory,  the  level  of  C'  is  increasingly  overpredicted  by  the 
transonic  method.  A  correction  for  the  effect  of  boundary  layer  and  wake 
for  section  3  at  M  =  0,82  (Pigs  8  and  9),  computed  with  the  LTRANV 

CO 

program,  partly  explains  the  overprediction  of  C'. 

z 

With  increasing  Mach  number  also  the  shape  of  the  spanwise  distribution 

of  C  is  predicted  with  decreasing  accuracy.  The  experimental  results 
z 

indicate  a  relatively  strong  influence  of  Mach  number  on  the  unsteady 
airloads  between  30  and  TO  %  of  the  semi-span.  The  theoretical  results, 
however,  show  a  typical  subsonic  distribution  of  which  only  the  level  is 
affected  by  increasing  Mach  number.  The  disagreement  compared  to  the 
transonic  experimental  results  may  indicate  shortcomings  of  the  subsonic 
3-d  correction  procedure  in  the  Q3D  method. 

Drawing  a  definite  conclusion  about  the  accuracy  of  the  Q3D  method, 
however,  is  difficult  because  of  the  inaccuracy  in  the  integration  of 
the  experimental  pressure  distributions,  as  mentioned  in  section  2.  This 
is  also  true  for  the  effect  of  amplitude  on  the  experimental  results, 
which  is  illustrated  in  figures  7  and  8  for  =  0.82.  Increasing  the 
amplitude  from  a 1  =  0.25  to  1  deg  leads  to  different  distributions  of 
quasi-steady  airloads  which  are,  probably  by  coincidence,  quite  similar  to 
the  results  of  subsonic  theory. 

The  sectional  moment  coefficients  (Fig.  9)  are  not  predicted  more 
accurately  in  a  quantitative  sense,  by  accounting  for  transonic  effects. 
However,  in  a  qualitative  sense  the  Q3D  method  correctly  predicts  the 
tendency  to  more  negative  (nose  up)  moment  coefficients  at  =  0.77, 
and  to  more  positive  (nose  down)  values  at  =  0.82.  This  tendency  is 
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explained  by  the  contribution  of  the  shock  wave  to  (V,  which  is  upstream 
of  the  quarter-chord  point  at  -  0.77  and  downstream  of  this  point 
at  M  =0.82  (Figs  6b  and  c).  Unlike  the  lift  coefficient,  the  moment 
coefficient  in  section  3  is  not  improved  by  applying  a  viscous  correction 
(Fig.  9,  =  0.82).  A  comparison  with  more  refined  calculation  methods 

is  necessary  to  determine  to  what  extent  the  small  perturbation 
assumption  and  the  qua si -3-d  concept  are  responsible  for  the  observed 
differences. 

Over»-i  1  coefficients :  The  above  observations  are,  of  course,  also  reflected 
in  the  overall  coefficients  in  figure  10.  It  should  be  noted  that  the 
overall  moment  is  not  an  average  of  the  sectional  moments,  as  it  is 
defined  relative  to  an  axis  through  the  quarter -chord  point  of  the  mean 
aerodynamic  chord,  normal  to  the  free  stream. 


5.2.2  Unsteady  results  for  f  =  2k  Hz 

Unsteady  results  at  a  mean  incidence  of  =  0.6  deg  and  a  frequency 
of  f  =  2k  Hz  are  shown  in  figures  11  to  ik. 

Pressure  distributions:  The  unsteady  pressure  distributions  at  f  =  2k  Hz 
(Fig,  11)  show  the  same  order  of  agreement  between  experimental  and 
theoretical  results  as  in  the  quasi-steady  case.  In  a  quantitative  sense, 
the  agreement  is  slightly  better.  The  computed  imaginary  parts  show 
clearly  the  increasing  influence  of  the  shock  wave  if  the  Mach  number  is 
increased. 

For  M  =0.82  the  computed  effect  of  amplitude  was  investigated. 

CD 

In  figure  11  the  significant  changes  in  pressure  are  indicated,  which 
occur  only  at  the  shock  wave  pressure  peaks.  The  Q3D  method  predicts  both 
a  reduction  and  a  slight  widening  of  the  peak.  The  experimental  results 
exhibit  only  a  reduction  of  the  pressure  peaks.  Here,  the  spacing  of  the 
pressure  holes  does  not  provide  sufficient  resolution  to  indicate  a 
widening  of  the  peak,  which  is  still  likely  to  be  present  in  reality  as 
in  the  quasi-steady  case.  As  a  result  the  computed  sectional  airloads, 
discussed  below,  are  not  significantly  affected  by  the  increase  of 
amplitude,  in  contrast  with  the  experimental  results  which  are  strongly 
dependent  on  such  a  change.  This  example  illustrates  very  distinctly  the 
remarks  in  section  2  on  the  experimental  accuracy  in  relation  with  the 
spacing  of  the  pressure  holes. 
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Sectional  coeff leients :  The  spanwise  distributions  of  unsteady  lift  and 

moment  coefficients  are  shown  in  figures  12  and  13,  respectively. 

The  lift  coefficient  behaves  very  much  like  the  quasi -steady 

case,  although  the  computed  effect  of  transonic  flow  is  smaller.  The 

agreement  between  the  theoretical  and  experimental  results  is  better  * 

in  spite  of  some  difference  in  the  theoretical  and  experimental  spanwise 

distributions  of  Ct. 

z 

The  computed  transonic  effects  on  CM  increase  with  Mach  number.  The 

z 

shift  of  C"  to  more  negative  values  is  correctly  predicted  by  the  Q3D 
method.  Although  the  experimental  values  are  somewhat  underpredicted, 
a  significant  improvement  has  been  obtained  compared  to  subsonic  theory. 

For  Mm  =  0.82  the  effect  of  an  increase  in  amplitude  to  =0.5  deg 
is  indicated  for  the  experimental  data;  in  the  computed  results  this 
effect  is  negligible. 

Also  the  distribution  of  C*  (Fig.  13)  shows  the  same  tendency  as  in 

the  quasi-steady  case  (Fig.  9).  The  computed  steady  flow  effects  yield  a 

more  positive  value  of  C"  at  transonic  conditions,  which  is  also  found  in 

m 

the  experimental  data.  In  spite  of  this  qualitative  improvement  relative 
to  linear  theory,  is  not  predicted  more  accurately  by  the  Q3D  method. 
The  overprediction  can  be  associated  with  the  overprediction  of  the 
unsteady  pressures  in  the  front  part  of  the  airfoil  {Fig.  11)* 

For  section  3  at  M  =  0.82  a  correction  for  viscous  effects  was 

oo 

computed  using  the  LTRANV  computer  program.  The  resulting  reduction  of 
(V  and  CM  (Figs  12,  13)  is  similar  too,  but  smaller  than  that  in  the 
quasi-steady  case.  The  effect  on  the  imaginary  part  is  a  positive  shift, 
which  causes  a  larger  difference  relative  to  the  experimental  results. 
This  again  points  to  shortcomings  of  the  inviscid  flow  computation, 
related  to  the  small  perturbation  assumption  and  the  quasi-3-d  concept. 
Overall  coefficients:  The  overall  unsteady  lift  and  moment  coefficients 
far  f  =  2U  Hz,  shown  in  figure  expectedly  reflect  the  behaviour  of  the 
sectional  airloads  discussed  so  far*  The  strong  variation  of  the 
experimental  unsteady  airloads  above  =  0.82  is  due  to  flow  separation, 
which  is  out  of  the  scope  of  this  report* 
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5.2.3  Unsteady  results  for  f  =  1+8  Hz 

The  differences  between  the  results  for  f  =  24  Hz  and  48  Hz  are 
caused  both  by  the  frequency  and  the  vibration  mode.  The  change  of  the 
vibration  mode  with  frequency  has  been  discussed  already  in  section  2, 

The  unsteady  results  for  48  Hz  are  shown  in  figures  15  to  18  for  a 
subsonic  (M  =  0.62)  and  a  transonic  (Mm  -  0.82)  case. 

Pressure  distributions:  The  unsteady  pressure  distributions  on  section  3 
(Fig.  15)  show  a  similar  development  with  Mach  number  as  already  discussed 
for  f  =  24  Hz. 

Sectional  coefficients:  Accounting  for  transonic  effects  leads  again  to 

an  increased  overprediction  of  and  to  a  better  prediction  of 

(Fig.  16).  The  maximum  values  of  the  experimental  lift  coefficients  at 

M  =  0.82.  which  occur  between  30  and  70  %  of  the  semi-span  like  in  the 

quasi-steady  case,  are  not  well  predicted. 

The  experimental  moment  coefficients  (Fig.  17)  are  predicted 

surprisingly  well  in  the  subsonic  case.  Apparently  the  errors  due  to  the 

neglect  of  viscous  effects  and  to  the  approximations  in  the  inviscid 

theory  cancel  out  here.  In  the  transonic  case,  the  Q3D  method  gives  a  good 

prediction  of  the  transonic  effects  on  (V  ,  which  are  not  predicted  by 

subsonic  theory.  On  the  other  hand,  the  Q3D  method  overpredicts  at 

M  =  0.82. 

00 

Overall  coefficients:  The  overall  unsteady  airloads  for  f  -  48  Hz  (Fig.  l8) 
show  a  behaviour  which  is  quite  similar  to  that  at  f  =  24  Hz. 


5.3  Unsteady  airloads  as  function  of  frequency 

From  the  results  for  a  =0.6  deg  presented  in  the  previous  sections 

m 

it  can  be  derived  that  the  effect  of  frequency  is  reasonably  well 
predicted  by  the  Q3D  method .  To  illustrate  this  effect  more  clea.ly , 
unsteady  results  will  also  be  presented  for  various  frequencies,  for  a 
mean  incidence  of  a  =2.6  deg,  at  M  =  0.77*  (Figs  19  to  22).  The  results 

IQ  ™ 

correspond  to  the  mean  steady  state  shown  in  figure  6d.  Again,  it  should 
be  noted  that  also  the  vibration  mode  varies  with  frequency. 

Pressure  distributions:  The  quasi-steady  and  unsteady  pressure  distribu¬ 
tions  on  section  3  shown  in  figure  19  exhibit  an  increasing  agreement 
between  theory  and  experiment  at  increasing  frequency.  In  particular 
the  decrease  of  the  real  part,  and  the  increase  of  the  imaginary  part 
of  the  pressure  peak  at  the  shock  wave  is  well  predicted. 
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Sectional  coefficients:  The  unsteady  airloads  are  shown  in  figures  20  and 

21.  In  particular  the  experimental  distributions  of  C'  and  C"  (Fig.  20) 

z  z 

show  strong  variations  over  the  wing  span,  which  are  likely  due  to  the 
experimental  inaccuracy  discussed  in  section  3.  The  nature  and  accuracy 
of  the  computed  distributions  is  rather  similar  to  that  of  the  previous 
cases,  and  need  not  be  discussed  further. 

Overall  coefficients:  The  development  of  the  overall  unsteady  airloads, 
shown  in  figure  22,  reflects  the  global  accuracy  of  the  Q3D  method  and 
demonstrates  that  the  merit  of  this  method  is  the  improved  accuracy  of 
the  computed  imaginary  part  of  the  unsteady  lift.  Note  the  discontinuity 
indicated  at  k  =  0.1,  which  corresponds  to  the  natural  frequency  of  the 
first  wing  bending  mode. 


Unsteady  airloads  as  function  of  mean  angle  of  attack 

The  effect  of  mean  angle  of  attack  is  investigated  for  =  0.77  at 

f  =  24  Hz.  The  steady  state  pressure  distributions  are  shown  in 

figures  6c  (a  =  0.6  deg)  and  6d  (a  =2.6  deg).  The  results  for  these 
o  o 

cases  are  already  discussed  separately  in  the  previous  sections. 

Pressure  distributions:  The  computed  pressure  distributions  (Fig,  23) 
predict  correctly  the  increase  and  downstream  displacement  of  the  shock 
wave  pressure  peaks  with  increasing  incidence. 

Sectional  coefficients:  The  computed  distributions  of  C'  (Figs  2k  and  25) 

show  an  increase*  while  Cf^  becomes  more  negative  at  larger  angle  of 

attack.  The  distributions  of  CT  and  C"  shift  slightly  to  larger  values, 

mm 

due  to  the  nose-down  contribution  of  the  more  downstream  positioned  shock 

wave*  The  above  effects  of  incidence  are  not  clearly  present  in  the 

experimental  results,  except  for  the  behaviour  of  C’  and  at  section  3- 

Overall  coefficients:  The  effect  of  incidence  on  the  overall  airloads  is 

shown  in  figure  26.  The  slight  effect  of  incidence  on  the  results  of 

subsonic  theory  is  due  to  a  small  variation  of  the  vibration  mode.  The 

increase  of  the  unsteady  airloads  with  increasing  incidence,  as  predicted 

by  the  Q3D  method,  only  agrees  with  the  experimental  results  for  CJ  and 

C"  In  a  quantitative  sense,  only  C"  is  predicted  satisfactorily,  as 
M  L 

before . 
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CONCLUSIONS 


Experimental  data  for  the  LANN  supercritical  wing  model  oscillating 
in  pitch  have  been  used  for  verification  of  a  quasi-3-d  method  to  compute 
unsteady  airloads  at  transonic  flow  conditions.  An  outline  has  been  given 
of  this  method,  which  combines  2-d  transonic  theory  with  2-d  and  3-d 
subsonic  theory*  A  comparison  of  computed  and  measured  results  for  a 
number  of  attached  flow  conditions,  covering  variations  of  Mach  number, 
frequency,  mean  angle  of  attack  leads  to  the  following  conclusions: 

a.  Accounting  for  the  effects  of  transonic  flow,  including  shock  wave 
motions,  using  the  quasi-3-d  method  does  not  in  all  respects  lead  to 
an  improved  prediction  of  experimental  unsteady  airloads,  compared  to 
subsonic  theory. 

b.  In  a  quantative  sense,  the  imaginary  part  of  the  sectional  and  overall 
lift  coefficients  is  predicted  significantly  better,  whereas  the  real 
part  is  more  overpredicted.  The  moment  coefficients  are  not  predicted 
more  accurately,  compared  to  subsonic  theory. 

c.  The  quasi-3-d  method  gives  a  qualitatively  better  prediction  of  the 
effects  of  Mach  number  and  frequency  on  the  unsteady  lift  coefficients 
and  on  the  real  part  of  the  moment  coefficients,  compared  to  subsonic 
theory.  The  predicted  effects  of  incidence  on  the  unsteady  airloads  do 
not  agree  very  well  with  the  measured  effects, 

d.  The  still  existing  differences  between  experimental  data  and  results  of 
the  quasi -3-d  method  can  be  partly  attributed  to  the  neglect  of 
viscous  effects,  to  the  small  perturbation  assumption  underlying  the 
2-d  transonic  flow  computations  and  to  the  subsonic  3-d  corrections. 

On  the  other  hand,  due  to  a  number  of  inoperative  pressure  holes  the 
experimental  lift  and  moment  coefficients  are  of  insufficient  accuracy 
to  draw  a  definite  conclusion  on  the  accuracy  of  the  quasi-3-d  method. 
For  this  reason,  also  the  measured  effect  of  amplitude  is  inconclusive. 

e.  The  computed  effect  of  amplitude  on  the  unsteady  airloads  on  a 
representative  wing  section  is  very  small. 
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A. 1  Definitions  and  sign  conventions  (taken  from  reference  1 ) 


Z,  Zi,  2  q 


M,  Mi,  M 


x 


c 


Note:  Coefficients  derived  from  the  zeroth  harmonic  component  of 
unsteady  pressure  signals,  are  indicated  simply  as  "steady" 
coefficients  instead  of  "mean"  coefficients.  For  a  pure  sinusoidal 
signal,  there  is  principally  no  difference  between  steady  and  mean 
coefficients . 

Steady 

Steady  pressure  coefficient: 


Steady  sectional  normal  force: 


Z  =  Cz  qm,  c;  Cz  =  -  /  (Cpm+  -Cpm_)d(^) 


o 


Steady  sectional  pitching  moment  about  quarter -chord 


M  =  Cm  q^.c2 ;  Cm  =  -  /  (Cpm+  -Cpm_)(^  -0.25 )d(^} 


o 


Steady  wing  normal  force : 


Z  . 
wing 


o  o 
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Steady  wing  pitching  moment  about  aerodynamic  center: 

Mving=  CM-q«‘CAC-S;  CH  =  '  s7~~  J  /  { °Vm+  “Cpm_ )  (— ^)d(|) c2  . dy 

AC  o  o 

Unsteady 

Unsteady  pressure  coefficient : 

F. 

Cpi  =  C'  +  i  C"  =  — - — 

P  P  Qco.ai 

Unsteady  sectional  normal  force: 


Zi  =  2tt  . q  .  %r.  Czi .  Ao.  ela)t 

^oo  2 


Czi  =  Cz  +  i  Cz  =  -  -  /  {Cpi  -Cpi  )d(-> 

TT  +  —  C 

O 

Unsteady  sectional  pitching  moment  about  quarter-chord: 


Mi  =  2ir.q  . (— )  .Cmi.Aa.e  ; 
00  ^ 


Cmi  =  Cm  +  i  Cm  =  -  -  f  (Cpi  -Cpi  )  (~  -0.25)d(-) 

7T  +  *  C  C 

O 

Unsteady  wing  normal  force: 

“wing  = 

CZi  =  CZ  +  1  CZ  ■  ■  is  II  (“PU  -Cpi.)a(|).c.dy 

o  o 

Unsteady  wing  pitching  moment  about  aerodynamic  center: 

c 

i j  *  0  /  AC^  S  .  iwt 

Mi,  =  2ir*q  *  (— — )  ,  Cu*  -  Aa.e  ; 

wing  ®  2  2  Mi 


1 


x-x 


Cw.=  Cf  +  i  CH.  -  “ 


'AC  o  o 


“  /  /  (Cpi  +  -Cpi_}(— ~^)d(^)  .c2  .dy 


Quasi -steady 

Quasi-steady  pressure  coefficient: 


Cpq  =  Cp 1  - 


q  Aa 
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Quasi-steady  sectional  normal  force: 


Zq  =  Czq.Aa  ; 


,  Cz{a  +Aa)  -  Cz(a  -Aa) 

Czq  =  C'  =  ±  ^ - 

n  Z  ir  2Aa 


Quasi-steady  sectional  pitching  moment  about  quarter-chord: 


Mq  ■-  2ir.q  .(— )  .Craq.Aa  ; 

00  C 


Cmq  =  C  *  — . 
’  ra  it 


_  Cm (a  +Aa)  -  Cm(a  -Aa) 
dm  jb 


2Aa 


Quasi-steady  wing  normal  force: 


% 


=  2ir.q  . C  .Aa  ; 
ing  “  2  Zq 

1  C  (a  +Aa)  -  C  (a  -Aa ) 
n  _  pi  _  2  ±  a  rc 

Zq  "  Z  "  ir*  2Aa 


Quasi -steady  wing  pitching  moment  about  aerodynamic  cent 


Mqw 


ing  2  '’2  MQ 

o  CM^am+Aa^  “  CM^am”A“^ 

c  =  C'  -  —  M  m _ M  m 

Mq  M  “  ti’  2Aa 
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Fig.  2  Wing  planform  and  measuring  points 
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Fig.  3  Test  conditions  for  LANN  model 
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Fig.  4  Influence  of  Mach  number  and  frequency  on  node  line  position 
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Fig.  5  Outline  of  quasi-3-d  method 
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Fig,  6  Mean  steady  pressure  distributions  on  section  3 
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Fig.  7  Quasi-steady  pressure  distributions  on  section  3  {a  =  0.6  deg) 
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Fig.  9  Quasi-steady  sectional  moment  coefficients  (ctQ=  0.6  deg) 
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Fig.  11  Unsteady  pressure  distributions  on  section  3  (a  =  0,6  deg,  f-  24  Hz) 
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Fig.  13  Unsteady  sectional  moment  coefficients  (a  =  0,6  deg,  f=  24  Hz> 
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Fig.  14  Unsteady  wing  lift  and  moment  coefficients  (ct^s  0.6  deg ,  f=  24  Hz) 
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Fig,  15  Unsteady  pressure  distributions  on  section  3  (<*0=  deg,  f  =  48  Hz) 
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Fig.  16  Unsteady  sectional  lift  coefficients  (ot  -  0,6  deg,  f-  48  Hz) 
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Fig.  17  Unsteady  sectional  moment  coefficients  (a  =  0.6  deg,  £=  48  Hz) 
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Fig.  18  Unsteady  wing  lift  and  moment  coefficients  (a  =  0.6  deg,  f=  48  Hz) 
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Fig.  19  Quasi-steady  and  unsteady  pressure  distributions  on  section  3  (a  =  2.6  degp  Mm=  0.77) 
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Fig*  21  Unsteady  sectional  moment  coefficients  (a  =  2*6  degp  0*77) 
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Fig.  22  Unsteady  wing  lift  and  moment  coefficients  {a  =  2.6  deg,  H  =  0.77) 
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Fig.  24  Unsteady  sectional  lift  coefficients  {M  =  0*77,  f=  24  Hz) 
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Fig.  25  Unsteady  sectional  moment  coefficients  (M^=  0.77,  f=  24  Hz) 
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Pig.  26  Unsteady  wing  lift  and  moment  coefficients  (Me=  0,77,  f=  24  Hz) 
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